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Abstract 



The quark delocalization and color screening model is used for a sys- 
tematic study of the effective potential between baryons in the u, d and 
s sector. The model is constrained by the properties of baryons and N- 
N scattering. The effective potentials for the N-N (IJ = 01,10,11,00) 
channels and the N-A and N-T, (I J = |l, |0, |l,|0) channels fit the 
N-N, N-A and N-T, scattering data reasonably well. This model pre- 
dicts: There are rather strong effective attractions between decuplet-baryons; 
the effective attractions between octet-baryons are weak or even repulsive; 
and the attractions between decuplet- and octet-baryons lie in between. 
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I. INTRODUCTION 



Realistic models for baryon-baryon interactions are relevant to the study of the strong 
interaction, the understanding of properties of nuclei and hypernuclei, multiquark states, 
neutron stars, stangelets, strange matter and so on. A simple but realistic baryon-baryon 
potential is the basic input for any realistic calculation of hypernuclear properties and the 
equation of state of strongly interacting matter . Such a potential is also needed for 
the study of dibaryons, as it is clear that an effective attraction between two baryons is 
necessary for the formation of a dibaryon, although it is not sufficient. 
There are several approaches to obtaining the effective potential between baryons. 
One-boson-exchange- models, such as the Nijmegen and Julich models [|J, (although there 
are many others,) are a typical approach. The advantage of this approach is that it gives a 
very accurate description of the nucleon-nucleon (NN) interaction. For extension to the 
hyperon-nucleon (YN) interaction, the problem of too many unkown parameters (coupling 
constants) can be solved by imposing SU3 flavor symmetry to constrain the coupling 
constants. The disadvantage is that a phenomenologically repulsive core, which is 
channel-dependent, has to be used. 

Another approach is based on applying the resonating group method (RGM) of 
calculation to quarks. Here, however, gluon exchange between the constituent quarks can 
only give the short range, repulsive core for the NN interaction. So in this model, meson 
exchange must be introduced to reproduce the medium and long range part of 
baryon-baryon interactions. It does have the advantage that the short range part is 
described by the quark model instead of a phenomenological repulsive core. 
Although the experimental data can put some constraints on models of baryon-baryon 
interactions, the existing data, even including the newest data, still can not discriminate 
among most models. This is true because any disagreements with the newer data can 
generally be repaired by fine-tuning parameters, relaxing some of the constraints or 
making other adjustments in the models ||. So in order to have a meaningful model for 
the baryon-baryon interaction with some predictive power, a model which is 
"parameter-free" is highly desirable. 

Recently a new quark model, the quark derealization, color screening model, (QDCSM) 
|J, has been developed. On comparison with the conventional quark cluster model, it is 
apparent that two new ingredients have been introduced. One is that the Hilbert space is 
enlarged to include fully confined and fully deconfined configurations as two extremes so 
that the actual configuration is determined variationally by the dynamics of the six quark 
system. The other is that the possible difference of the q — q interaction inside a hadron 
(in an overall colorless state) and between quarks from two different colorless hadrons, due 
to the nonlinearity of QCDQ, is taken into consideration. 

In this way, without introducing mesons, the long standing problem of the missing medium 



For example, the contribution from gluons exchanged among three quarks which are in a color- 
singlet state is zero, while there is a nonzero contribution among three quarks which are in "colorful" 
state. 
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range attraction of NN forces in quark models is solved, and the parallel between nuclear 
and molecular forces obtains a natural explanation. The model has been successfully 
applied to describe NN scattering PJlQfl . Using this model, a systematic search for dibaryon 
candidates within the three flavor world of u, d, s quarks was also made [|TT | . A relativistic 



model calculation obtains qualitatively similar results, even though the confinement 
mechanism differs. Therefore, it seems worthwhile to apply this model to perform 
systematic and parameter free calculations of the effective baryon-baryon potentials. 
This paper is organized as follows: In section II, the model Hamiltonian and Hilbert space 
are described. Section III is devoted to a sketch of the calculation method. The results are 
given in section IV and a conclusion is given in section V. 



II. QUARK DELOCALIZATION, COLOR SCREENING MODEL 

The details of the model can be obtained from ref. |T^,|n|]. Here we only give a brief review. 
For a single hadron, the QDCSM is same as the quark potential model, which has been 
successfully applied to describe the properties of baryons and mesons. The model 
Hamiltonian which describes a single baryon is written as 



H(3) = ^(m, + 



i=i 



Pt 
2m,- 



V- ■ = V c + V G 



V; 
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The rrii, pi and r»j are quark masses, momenta and pairwise separations. The Aj are color 
SU(3) Gell-Mann matrices and a s is the strong coupling constant. For the confinement 
potential, V c , both linear and quadratic forms are used in our calculations. In the effective 
one gluon exchange potential, V G , only the color Coulomb and color magnetic terms are 
retained; other terms have been neglected temporarily |10| in order to reduce the 
calculational burden. 

The wave function (WF) for a single baryon has the form 



^(123) = X (123)r /5/J (123)0(123). 

Here x(123) is the singlet color WF, ^5/j(123) is the symmetric spin-flavor 
SU2X/ DSU/xSUg WF (with S =strangeness, I =isospin, J =spin). The spatial 
wavef unction is given by 

0(123) = 0i s (ri)0i s (f 2 )0i s (f 3 ), 



1 \ 4 (r-s) 2 



nb 2 



e a*- 3 



(2) 



(3) 



where s is a reference center and b is a baryon size parameter. 
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For the two baryon system, we extend the Hamiltonian to 



H(Q) = J2(m t + A.) + £ - T c (6), (4) 

i=l 1 i<j=l 

but modify the color confinement as follows: 

^ .,. ■ —\i ■ XjCtr^j if i, j occur in the same baryon orbit, 

■ — A» • A, — (1 — e~^ nT ^) if i.j occur in different baryon orbit. 

Screened linear color confinement is found in lattice QCD calculations JHJ and can be 
parametrized as 



V lr) = + —j, (6) 

with a s = 0.21 ±0.01, v^ = 400MeV, ^ l = 0.90 ± 0.20 fm. 

For the two baryon wave function, we extend the quark cluster model space by introducing 
delocalized single quark orbits: 

<Mr) = (<Mr) + e(s)(f> R {f)) /N(s), 

= {<Pr(t) + e(s)0i(r)) /N(s), (7) 
N 2 (s) = l + e 2 (s) + 2e(s)(<p L \<f )B ). 

where 4>l,<Pr are the quark cluster bases of the type described above for individual baryons: 



^Ly) — ( ^2 J e 2^ (left centered orbit), 

^(rj — (^2 ) e ^ (right centered orbit) 



Here s is the separation between the centers of two q 3 clusters. The derealization 
parameter e(s) is determined variationally for every s = \s\ by the q 6 dynamics (see section 
III). The quark molecular orbit introduced by Fl. Stancu and L. Wilets |TJ[] is a restricted 
version of this structure. 
The two baryon wave function is written as 

K lFl , a2 F 2 a • • "6) = ^[Vw^(123)VW456)] a , 

^(123) = x(123)r /5l/lJlFl (123)^(l)^(2)^ L (3) (9) 
(456) = X (456)r/ 52/2j2F2 (456)^(4)^(5)^(6). 

Here a = {SIJ) are again the strong interaction conserved quantum numbers: strangeness, 
isospin and spin. (Orbital angular momentum is assumed to be zero for the lowest states. 
In principle, an angular momentum projection should be done but we leave this for future 
refinement.) The q 3 cluster WF is almost the same as given in Eq.(0), but the single 
cluster Gaussian WF, Eq.(^j), is replaced by the delocalized orbital WF, eq.(^). Finally, A 
is the normalized antisymmetry operator 

A = -^V(-)V 
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The model parameters m, m s , b, a s , a are fixed by the spectrum of baryons; fix is taken 
from lattice gauge QCD calculations; and fi% is fixed by fitting the N-N phase shifts. The 
fitted parameters are given in Table I. 



Table I. The parameters used in the calculations. 





m(MeV) 


m s (MeV) 


6(fm) 


a s 


a(MeV-fm- 2 ) 


H"n 


n = 1 


313 


618.45 


0.625 


1.71 


39.14 


1.1 


n = 2 


313 


633.76 


0.603 


1.54 


25.13 


1.6 



III. CALCULATION METHOD 

Since the relative motion between two baryons is rather slow compared to the internal 
motion of the quarks in the baryons and the quantum fluctuation of the center of mass 
coordinates of the baryons is neglected temporarily (again leaving refinement for the 
future), we can define the effective baryon-baryon potential as follows: 

V a (s ) = E 6 (s = s ) - E e (s = oo) (10) 

where Eq(s) is the energy of two baryon system at separation s, which can be obtained by 
calculating the matrix elements of the Hamiltonian, Eq.(Q), in the baryon states, Eq.(|9]). 
In the nonrelativistic approximation, 

E 6 {s) = £ 3 (123) + £ 3 (456) + E rel + V a {s), (11) 

where E3 is the total energy of a single baryon and E re [ is the energy of relative motion of 
two baryons. When s — > 00, we must have V a — > 0, so 

E 6 (s = 00) = £73(123) + £ 3 (456) + E rel . (12) 

Note that the six-quark WF decomposition maintains a non-zero E re i even at 00 and it is 
our assumption that this constant is the appropriate quantity to remove at finite s which 
leads to the expression for V a (s) shown in Eq . (|T0l) . 

The treatment of the center of mass motion merits some additional discussion. In order to 
have a correct separation of total kinetic energy into internal, relative and center of mass 
parts in the asymptotic region, we redefined the CM energy using an average of the single 
quark mass 

TM = ^ (±4 = i (^i + ±±) (±P,)\ (13) 
2M n J 2 V n m n mj J 

where n and n s are total number of quarks and s-quarks in the system, respectively. In 
this way, the kinetic energy is free from CM motion when an SU(3)-flavor symmetric wave 
function is used. 

T(n) = ±^P1- T c {n) = J- (V^L + ±™) ± (p { - p 3 f (14) 
~[ 2rrii 2nm \ n n m s J i^J=i 
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A similar treatment can be found in ref. 0. 

To calculate the matrix elements of H in the basis given by Eq.(||) is a daunting task, 
especially for a systematic search, where many matrix elements must be evaluated. So a 
group theory method was developed to relate the physical basis to a group chain 
classification basis (symmetry basis). In the symmetry basis, the well known fractional 
parentage expansion can be used. The calculation of six-body matrix elements reduces to 



two-body matrix elements and four-body overlaps. For details of the method, see ref. [T5 
In order to discern the effect of channel coupling, a channel potential is introduced. For a 
given set of quantum numbers a = (SIJ), the effective channel potentials are defined by 
the same formula Eq. (|T0|) , where E 6 (s = s ) and E$(s = oo) are the eigenenergies of the 
two baryon system at separation so and oo, respectively, which can be obtained by 
diagonalizing the q 6 Hamiltonian. The eigenstates are expressed as multiple physical 
channel coupling wave functions 

*„(1...6)= £ C: iFua2F2 K lFua2 F 2 - (15) 

ctiFi,at2F2 

The channel coupling coefficients C% F a2 p 2 are determined by the diagonalization. 
In general, we can identify which of the baryon channels an eigenenergy corresponds to by 
finding which baryon channel for that eigenenergy has the largest channel coupling 
coefficient. Although this identification is not unique, it guarantees that the channel 
energy varies smoothly with s - Where no coefficient is discernably dominant, we examine 
the energies at (nearby) larger and smaller values of so and again choose the value which 
provides for a smooth continuation of the channel energy. This replaces an earlier 
procedure we used which involved choosing the lowest eigenenergy available. This older 
procedure had the disadvantage of not having a clear relation to particular baryon-baryon 
channels, of introducing sharp (discontinuous in the derivative with respect to Sq) changes 
in the energy, and of prohibiting the identification of any but the lowest energy channel. 
The energies obtained in this way depend on the separation s and the derealization 
parameter e(s). We repeat the calculation for each s by varying e(s) from 0.01-1.0 with 
0.01 stepsize to obtain the minimum of the lowest eigenenergy. At the same time, the 
derealization parameter e(s) and all eigenenergies are determined. 



IV. RESULTS 



The effective potential between combinations of octet- and decuplet-baryons (168 channels, 
altogether) were calculated. To save space, only a few typical examples are shown in 
Figs. 1-7. (We will be happy to provide any of the others upon request.) 
Fig. la shows the effective N-N potential for channels (IJ) = (01), (10), (11) and (00) with 
the quadratic color confinement interaction. Clearly there is an effective attraction in the 
first two channels, with (01) a little stronger than (10). The other two channels are 
repulsive. These potentials are quite similar to phenomenological N-N potentials such as 
the Reid soft core N-N potential. [|16j However, when we proceed to a dynamical 
calculation of the phase shifts, we must also perform a partial wave decomposition and 
make a center of mass correction of the kernel of the interaction Hamiltonian. The net 
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effect of these calculations is to create some additional attraction, so that a slight 
reduction of the value of the screening parameter is then needed to produce good 
agreement with the results of N-N phase shift analyses. |T7 



Fig. lb is similar to Fig. la but with the linear color confinement interaction. Comparison 
of Fig. la with Fig. lb makes it clear that both sets of potentials have qualitatively the 
same pattern, but with a weaker attraction for the case of linear confinement. In fact, we 
find that for the linear color confinement interaction, the potentials in these four channels 
agree reasonably well [D| with the phase shift analyses without any change from the lattice 
QCD value of the screening parameter JOJ. In this sense, our results provide a 
parameter-free prediction of the phase shifts in these channels. 

The channel coupling calculation shows that channel coupling produces minor effects for 
the N-N channels, and so we do not present those results here. In the following, we only 
present the results with quadratic color confinement to save space. 

Figs. 2-3 shows the results for AT-E and N-A channels. For N-A, the spin triplet state is a 
little more attractive than the spin singlet; the channel coupling adds a bit more attraction 
to this state but leaves the spin singlet almost unchanged. We note that and \He both 
have spin zero ground states and spin one excited states. One might interpret this as 
evidence that the spin-singlet iV-A interaction is more attractive than the spin triplet. 
However, the situation is not so transparent due to the complications of the interactions of 
the four bodies involved and, in addition, A-E° mixing effects. The spin one ground state 
of the deuteron is certainly an indication that the spin triplet N-N interaction is more 
attractive than the spin triplet, and we might reasonably expect (by flavor symmetry) that 
this should hold true for all octet-baryon combinations. However, due to the strong tensor 
interaction from pion exchange, the Nijmegen OBE model F [|J nonetheless includes a 
more attractive spin singlet N-N central interaction. Furthermore, there is a paucity of 
direct data on scattering in the Y-N channels, and widely differing relative strengths for 
the central interaction are all consistent with both the available data and the nuclear 
states referred to above. The QDCSM, on the other hand, predicts that spin triplet N-N 
and A^-A interactions are stronger than spin singlet ones. Clearly, which central interaction 
is stronger in each case merits additional study. 

For A-E, we find the strongest attraction in the IJ = |l channel, while the I J = |0, and 
|0 channels both have a little weaker attraction (single channel case), and the I J = |l 
channel is repulsive. Channel coupling has little effect on ~1, and pushes |0 from 
attractive to repulsive. These show that the Af-E potential is more strongly spin and/or 
isospin dependent than the A^-A potentials, which have a little weaker dependence on spin. 
These results are in qualtitative agreement with the calculations of OBE models M and 
hybrid quark model calculations 0, except that our attraction for the A r -E(|l) channel is 
too strong. It is worth noting that the Af-E (|l) channel has a small bump at R s ~ 1.3fm. 
The reason for this is not clear, but it seems to correspond to the statement of ref. ||, that 
the repulsion between two clusters has its maximum not at distance zero, but at a finite 
cluster distance. 

For all other channels, some general features of the effective potential have been found in 
this calculation: An effective attraction exists in most channels. The attraction is very 
strong between two decuplet-baryons, but the attraction is very weak and even repulsive 
between octet-baryons, while the strength of the attraction between octet- and 
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decuplet-baryons lies in between. The reason is that in a decuplet baryon, spin is totally 
symmetric and color is totally antisymmetric, so that the hyperfine interaction is purely 
repulsive. For a six quark system, the spin and color are neither totally symmetric nor 
antisymmetric, leading to the hyperfine interaction not always being repulsive, and so 
providing an effective attraction. 

Some states are forbidden for structureless baryons by the Pauli principle: A A (/J = 33), 
A£*(§3), £*H*(§3), fifi(03). For these states, we find that the energy of the system 
becomes very large when the distance between the two baryons tends to zero (see Fig.4). 
Because the effective B-B interactions are directly related to multiquark states, in the 
following we present the baryon-baryon potential along with a discussion regarding 
dibaryon states. 

(a) States with only non-strange quarks. Since the states in this part do not contain any 
strange quarks, they are not affected by symmetry breaking. 

Clearly many channels have an effective attraction, but to form bound states one has to 
consider not only two body decay modes, but also many body decay modes. For example, 
for a AA channel to form a narrow dibaryon resonance, the energy of dibaryon should not 
only be less than the sum of the masses of two As, but it should also be less than the sum 
of the masses of two nucleons and two pions because each A itself can decay into Nit. In 
addition, the zero-point harmonic oscillater energy must be added to the system. Only 
channels with sufficient attraction can form bound states or narrow resonances. 
The interesting states are: 

NN (01), NN (10) and AA(03). The NN (01) state is a loosely bound state: This is just 
the deuteron, which is reproduced here. The NN (10) state corresponds to the known NN 
zero energy resonance. The A A (03), called d*, is a narrow resonance. Although the 
AA(01) channel is more attractive than the AA(03) (see Fig. 5), it couples strongly to 
NN (01), so it is not a good candidate for an observable dibaryon. Conversely, because the 
A A (03) has a large angular momentum, it cannot couple to an NN s-wave state, or to 
NA (also forbidden by isospin) directly. The large angular momentum hinders the decay 
to NN, and the energy of the state prevents its decay to NNtttt, so it is a good candidate 
for a narrow resonance. |2D| 

(b) States with one or more strange quarks. The general features stay the same for states 
with strange quarks. There are also several particularly interesting states (see Figs. 6-8). 
££(20), SS(00) are both deuteron-like states, consisting of two octet-baryons, which have 
a weak attraction and the minima of their potentials are found at rather large separations 
(> 1 fm). Channel coupling has a minor effect on the results. These are possible dibaryon 
candidates but are sensitive to model details. 

It is worth mentioning the fif2(00) state. In our calculations, its effective attraction is not 
weak, but rather strong enough to form a strong interaction bound state. Because it 
consists of two Qs and Q is strong interaction stable, it is a good candidate for a dibaryon, 
although it is also sensitive to model details. In addition, it may be too difficult to produce 
this state even in relativistic heavy ion collision experiments. 

The H particle (YIJ = 000) is a special state in our calculation. The channel coupling 
calculation (see Fig. 6) shows that it is a six quark state consisting mainly of octet-baryons 
(dominant components are iVS and ££) that has an attraction that is not too weak. It is 
possible that it appears as a dibaryon. 
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The states YIJ = -1±2 and both have a rather strong attraction. However, when 

channel coupling is taken into account, the attraction is greatly reduced for the J = 1 
channel, whereas the J = 2 channel is almost unaffected. Although they consist of octet- 
and decuplet-baryons, their energies are lower than the AEn three body decay channels. A 
dibaryon may form in the state YIJ = — 1|2. 



V. CONCLUSION 



Taking into account both color screening and quark derealization, a systematic study of 
the effective baryon-baryon potentials has been performed by using a quark cluster model. 
The model parameters were fixed by baryon properties and the color screening parameter 
was taken from lattice gauge QCD calculations for linear confinement (or from N-N 
scattering data for quadratic confinement). Consequently, the calculated baryon-baryon 
potentials can be viewed as a parameter free result. As a general trend, we find that a very 
strong attraction exists between two decuplet-baryons. For decuplet and octet baryons, 
one has a mild attraction. Very weak attraction, or even repulsion appears in the 
octet-octet baryon system. 

Because the "predicted" N-N, N-A and N-H effective potentials fit the scattering data 
reasonably well, we expect that the other baryon-baryon effective potentials are reasonable 
predictions as well. At least for those cases for which one has neither any experimental 
data nor good theoretical models, our model results should be useful intermediate guides 
for studies of such exotics as mutliquark states, strangelets and of the equation of state of 
strongly interacting matter in relativistic heavy ion collisions, where multistrange baryons 
constitute a significant component. 

For a dibaryon search, these results suggest that experiments should focus on octet-octet 
systems with low spin and decuplet-decuplet systems with high spin. For octet-octet 
systems, there are dibaryon candidates which are loosely bound but stable against strong 
interaction decay and have small quark derealization. This is the nuclear type of dibaryon 
and the deuteron is a typical example. For decuplet-decuplet systems, the most striking 
resonances are states with high spin and large quark derealization. These are a quark 
matter type of dibaryon and the d* is a typical example. 

Strangeness conserving transition potentials (e.g., iVA — > NT,) can also be obtained in our 
model. We defer discussion of those results to a later paper. 

We thank R. Timmermans for a useful conversation. This research is supported in part by 
the Department of Energy under contract W-7405-ENG-36 and in part by the National 
Science Foundation of China. 
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FIGURE CAPTIONS 

Fig. la Effective potential in MeV vs. baryon separation in fm for N-N channels with 

quadratic confinement. 
Fig. lb Same as in Fig. la, but with linear confinement. 

Fig. 2 Effective potential in MeV vs. baryon separation in fm for iV-E channels with 

and without channel coupling. 
Fig. 3 Effective potential in MeV vs. baryon separation in fm for N-A channels with 

and without channel coupling. 
Fig. 4 Effective potential in MeV vs. baryon separation in fm for various forbidden 
states. 

Fig. 5 Effective potential in MeV vs. baryon separation in fm for A-A channels. 

Fig. 6 Effective potential in MeV vs. baryon separation in fm for YIJ = 000 channels 

(i.e., with two strange quarks) with and without channel coupling. 
Fig. 7 Effective potential in MeV vs. baryon separation in fm for N-Q, channels with 

and without channel coupling. 
Fig. 8 Effective potential in MeV vs. baryon separation in fm for EE, SS and flfl 

channels. 
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Fig. 1a. Effective potential 

(for NN with quadratic confinement) 
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Fig. 1b. Effective potential 

(for NN with linear confinement) 



X 

\ 




I 

IJ=01 


\ 




IJ=10 


\ 




IJ=11 


~ \ 

\ \ 

X 

\ 


\ 


IJ=00 


\ 

\ 


\ 




\ 

\ 


\ 




\ 

\ 


\ 




X x 

\ 


\ 




\ 


\ \ 


- 


\ 

\ 

\ 


\ \ 




_ \ 

\ 


\ \ 

\ 




^^^^ \ 

\ 


\ \ 






\ 

\ \ 

\ \ 






\ \ 

\ \ \ 

\ X ^ 






\ X "\ 

\ X 

\v X 

\ \ 




I 



s(fm) 



Fig. 3. Effective potential 

(for NA) 
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Fig. 4. Effective potential 

(forbidden states) 
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Fig. 5. Effective potential 
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Fig. 6. Effective potential 

(for YIJ=000) 








s(fm) 



Fig. 7. Effective potential 

(for HQ) 
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Fig. 8. Effective potential 

(for EE, EE, 
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